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O . Summary 



Berry and Klein (J. Mod. Opt. (1997) 43 2139-2164) showed that the Talbot effects 
in classical optics are naturally expressed by Gauss sums in number theory. Their result 
was obtained by a computation of Helmholtz equation. In this article, we calculate the 
^ ' effects using Fresnel integral and show that the result is also represented by Gauss sums. 
^ However function forms of these two computational results are apparently different. We 
^ show that the reciprocity law of Gauss sums connects these results and both completely 
. agree with. The Helmholtz equation can be regarded as an equation based upon wavy 
O I nature in optics whereas the Fresnel integral is defined by a sum over the paths based 
^ upon a particle picture in optics. Thus the agreement of these two computational results 
^ could be interpreted in terms of the concept of the wave-particle complementarity, though 
the concept is for quantum mechanical phenomenon. This interpretation leads us to a 
fH [ relation between the reciprocity of Gauss sums in number theory and the wave-particle 



complementarity in wave physics. We discuss it in detail. 
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§1. Introduction 

The Talbot effects are diffraction grating phenomenon in classical optics, discovered by 
Talbot (1838), which are known as self-image (Patorski 1989); due to the effects, same 
patterns arrayed on one or two dimensional grating sharply recover on the screen for a 
certain condition without lens system. It is remarkable that such the effects occurs due to 
wavy properties. 

Recently Berry and Klein (1996) discovered that behind the phenomenon there is an 
arithmetic structure. They considered optical wave with wavelength A which stems from 
a one-dimensional grating with a period a. The grating is set up so that its plane agrees 
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with the wave front and the screen is set parallel to the grating with distance z. They 
investigated Talbot effects when z was a fractional number times the Talbot distance 
Zt '■= 0? /\. They found a beautiful connection between wave physics and number theory 
in the investigation of the Talbot effects, which is based upon a deep insight by Hannay 
and Berry (1980). Berry and Klein investigated the Talbot effects by solving the Helmholtz 
equation (Noponen k, Turunen 1993) and showed that the Talbot effects were explained 
by Gauss sums and the quadratic reciprocity law (see for example Ireland & Rosen 1990, 
Chap. 6). 

In this article, we present another connection of the problem with Gauss sums based 
upon the arguments of Berry and Klein. We will also neglect the polarization effect (Nopo- 
nen & Turunen 1993) and partially review investigation of Berry and Klein by solving the 
Helmholtz equation in §2. In §3, we compute the same system by means of Fresnel integral 
(Winthrop & Worthington 1965), whose apparent function form differs from the results in 
§2. However in terms of the reciprocity of Gauss sums we will show the agreement of both 
results in §4. We should note that whereas Helmholtz equation is a wave equation, the 
Fresnel integral is an integration over the optical paths with a natural integration weight. 
The former one should be regarded as an expression of the wavy properties and later one 
should be interpreted as a transformation to wave expression from the particle nature. 
Thus the agreement of the both results reminds us of wave-particle complementarity, i.e., 
the wave-particle duality of nature in the quantum mechanics (Bohr 1928). We also argue 
that the complementarity is related to the reciprocity of Gauss sums (Hecke 1981 Chap. 8) 
in the Talbot effects, by using the analogy between the quantum mechanics and the optics 
in §4. 

In order to make the agreements between two computations from Helmholtz equation 
and Fresnel integral easier, we sometimes employ different expressions from those in the 
papers of Berry and Klein and Hannay and Berry (loc.cit.). Thus the difference between 
our formula and theirs sometimes occurs in this paper. 

Interestingly in the same period as Talbot's discovery. Gauss lived and studied optics, 
number theory including Gauss sums and quadratic residues and so on (Klein 1926). 



§2. Wavy Expressions 

As in the paper of Berry and Klein (1996), we will start with an incident plane wave of 
wavelength A coming through periodic S functions transparency whose period is a. Such a 
grating is called S-comb. Let the transverse direction denote x while z denotes along the 
optical axis. 

In this section, we will deal with the system in terms of the Helmholtz wave equation 
following Berry and Klein (1996) and Noponen and Turunen(1993). Let us introduce the 
dimensionless parameters with respect to Talbot distance zt, 



On the 5-comb grating plane ( = 0, the optical wave is expressed by the Poisson sum 
formula, 

oo oo 

^com6(e,0)= Yl exp(27rien)= Yl ^i^'^)- (2-2) 

n=— oo m=— oo 

As there is a discrete translation symmetry, we can set ^ e R as ^ = Co + n for £ 
(—1/2, 1/2) and an integer n. By letting 

oo 

il^comb{C,0= Yl ^(C)exp(27ri^n), (2-3) 

n= — oo 

we substitute it into the Helmholtz equation, 

+ ^+\^] I V'comb = 0. (2-4) 



Then we obtain the solution of (2-4) up to a constant factor C, 



oo / 

i^combiC, = C Y exp(27rzCn) exp j 27rzC (^^) 



1 - 



nX 



a 



(2-5) 



Here we note 5{cx) = 5{x)/c. For the case n > a/ A, the exponential factor becomes a 
damping function. 

By employing the paraxial approximation, (2-5) is expressed by. 



-ipcombiC, C)^C Y exp(27ri^n) exp ( 27riC (^^) 



(2-6) 



Of course, large n is out of the approximation but as it is expected that the contribution 
from the large n is less than that from small n in Fourier analysis, we will go on to consider 
(2-6) as in the paper of Berry and Klein (1996). Let the right hand side of (2-6) be denoted 
by il!p{C, () exp{ikz) following the notations in the paper, where k := 27r/A. We obtain, 

MC, = C Y exp(27rzen) exp (^-277^ (^nA ) . (2-7) 

n=— oo ^ \ / / 

Let us consider the case that C is a rational number, i.e., for coprime positive numbers 
p and q 

C = (2-8) 

By introducing the integers I, s such that n = Iq + s, the quantity of each term in (2-7) 
become, 

Cn - = ^s-^^s^ + ilq - ^l^qp - Ips. (2-9) 

3 



Further we note the relation, 



^-■niqpl _ ( — 1)5' ^ = e" 



niqpl 



and it is the unity if pq is even. Then we have 

q-l 



q — l oo 

s-ps^)/q g7ri(2Cg-pg)Z 



s=0 l= — oo 

Using the Poisson sum formula again, we obtain. 



oo ^ 



(2-10) 



(2-11) 



(2-12) 



where 



q-l 



1, if odd, 
0, if even. 



1 

A{n] q,p) = exp {in [{2n + qcqp) s - ps^] /q) 

v5 s=0 



(2-13) 
(2-14) 



We note that if pq is odd, {2n + qcgp) is also odd. Hence if we can use the formula given 
by (A-18) in the Appendix which is based on the Appendix of the paper of Hannay and 
Berry (1980), then (2-14) is expressed as follows (Berry & Klein 1996); 



exp ZTT 









"1" 










) 


.p. 







A{n;q,p) = < 



P 



exp tn 



exp — «7r 



1 p 
4 q 



n 



4 q 



2p 



(2n + qf 



where 



is a unique positive integer smaller than q satisfying 



P 



1 mod q, 



p even, q odd, 

p odd, q even, 

p odd, q odd. 
(2-15) 

(2-16) 



and is the Jacobi symbol which is a product of the Legendre symbol for the 



prime factors .s of q (Ireland & Rosen 1990, Chap. 5), 



+1, if there is an integer m such that = p mod s, 
— 1, otherwise. 
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(2-17) 



§3. Particle Expressions 



In this section, we will consider the system in terms of the Presnel integration following 
the study by Winthrop and Worthington (1965). Here we note that the Fresnel integration 
should be considered as a simple path integration (Schulman 1981, Chap. 20; Sanchez & 
Wolf 1985) over the shortest optical paths obeying the minimal principle of Fermat and 
these paths should be regarded as a particle picture in optics. 

Using the Fresnel integral for 5-comb wave function (2-2), we obtain the if^comb the 
point (^,C), 

i'comb{C,0 = J _ ^^t'^, ^ exp (j^Vi^ - ^ra' + z^^ i^combie, 0) (cos^ + 1), 

(3-1) 



where cos e = z/y/{C - O^a'^ + z'^. Here we note that (3-1) is a solution of (2-4) in the 
sense of (9 (A/ a/(^ — ^')'^a^ + z'^) as the function is an approximate kernel of the Helmholtz 
differential operator, called Fresnel approximation [BW, Chap. 8]. 

The optical distance between a point (^, Q and the n-th split (— n, 0) is given by 

^V(^ + n)%2 + ^2^ (3.2) 
and in the paraxial approximation, it is approximated by 

In this approximation, the denominator of the first factor is approximated by z and cos 6 
can be regarded as the unity, (3-1) becomes 

V'com6(e,C)~C" exp(27ri^/A)expf7rd ^ + + j , (3-4) 

n=— 00 ^ \ ' ^ S / / 



using a certain constant factor C . Of course, this approximation is not contradict with 
the Fresnel approximation . By letting the right hand side of (3-4) denoted by '0p ex.^{ikz), 
we have 



This formula is essentially obtained by Winthrop and Worthington (1965). In the formula, 
we have also assumed that contribution from large n is less than that from small n and 
have applied the paraxial approximation to large n case. (This application can be justified 

by the considerations that for the case that either ^ or ^ is an irrational number, the phase 
at large n becomes random phase for each n and sum of each terms with the random phase 
damps its amplitude and for the case that both Q and ^ are rational numbers, the net effect 
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from (3-1), in which the larger n contributes less, repeats to appear in (3-5) but rescaling 
C like (A-1) in the Appendix, (3-5) brings us the data of the net effect.) 

We also consider the same situation as (2-8) or ( = p/q for coprime positive numbers p 
and q. By introducing the integers I, s such that n = lp + s, the quantity of each terms in 
(3-5) becomes. 



{Qn+ -n t,s- + --S + ^Iq + -I qp + Iqs. 



2 p p 2p 
Noting (2-10), (3-5) is expressed by 

q-l 



s=Q 



Corresponding to (2-12), we obtain the wave function of the system. 



(3-6) 



(3-7) 



4(^>C) = C'VP £ A{n;q,p)5{i-\eqp--), 

A q 



n= — oo 



(3-8) 



where 



1 

A{n;q,p) := —^^^^P (^tt [{2n + qcqp) s + qs"^] /p+ (2n + qcqpf / Apqj . (3-9) 



As (2-15), we have 



p 



exp it: 



4^ 



V 



qp 



A{n;q,p) = < 



p 



CXp —ITT 



4 I p 



n 



qp 



n 



exp —m 



4 \ p 



2q 



Aqp 



{2n + qy 



p even, q odd, 

p odd, even, 
p odd, q odd. 



(3-10) 

We note that the function form in (3-10) differs from that in (2-15). Especially the apparent 
form of the final case of odd pq is completely different from that in (2-15); the roles of q 
and p in both formulae look inverted. 



§4. Discussion 



As we have two expressions of the same Talbot effect, we will discuss both expressions. 
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Due to Hecke (Hecke 1981 Chap. 8), we have the reciprocity of Gauss sums, if ab even 
and c is even (or if ab is odd and c is odd) (Hannay & Berry 1980), 



o-l 



at 
J 



1/2 



6-1 



, c X ^ 

\ Q ~F 



(at'^+ct) 



(4-1) 



t=o 



Here we use the convention \/i — i^l'^ = e^'^^^. Prom the function forms in (2-15) and 
(3-10), we have the relation. 



A{n; q,p) = ViA{n;q,p). 



(4-2) 



In other words, the result which we used the Helmholtz equation and that used by the 
Fresnel integration agree with after applying the paraxial approximation. Whereas we 
used the Fourier transformation on solving the Helmholtz equation, the computation of 
Fresnel integration is directly connected with the optical paths. As (4-1) is obtained by 
the Fourier transformation as mentioned in the Appendix (A-15) following Appendix in 
the paper of Hannay and Berry (1980), the correspondence is very natural. However as 
in both approaches, we used the essentially same but different paraxial approximations, 
(4-2) is not trivial. Further by considering the meanings of (4-2) as follows, it will turn 
out non-trivial. 

Hereafter let us consider the relation (4-2) primitively in order to reveal the arithmetic 
structure of (4-2). First we consider n-dependence up to constant factor. Since p and q 
are coprime. 



"1" 




"1" 






q + 




p 


.1. 


p 


.p. 





p ^1 + pq, 



and thus 



q + 



p ] = 1 + {pqY mod 2pq. 



(4-3) 



(4-4) 



Hence the correspondence of the n-dependence parts of even pq case is expressed by 

2 \ / \ 2 




2 _ P 

n = - 



n 



mod 2. 



(4-5) 



Similarly, for the case of pq odd, we obtain, 

2 \ 

(2n + 9)2 




{2n + qf + {^{p + q + l) + ^] {2n + qf mod 2 



(4-6) 



The derivation of (4-6) needs heavy computations in the sense of primitive number theory 
due to the factor 1 /4. We used the relations for pq odd case, 



"1" 
2 



1+p 





"1" 




"1" 




K 




p + 








p 


.P. 





(1 +pq) = ^ 1^ ^ mod 2, 



(4-7) 



which is proved by by expressing p — 4^ + a and q — + (3 {/j,, u are positive integers 
and a and /? are 1, or 3.) and by considering the possible cases. 

With respect to the n-dependence, except constant phase and factor, we have the same 
dependence. It implies that they express the same intensity of light, square of wave func- 
tion, in the screen. 

Next we will consider the equality of (4-2) including the constant phase factor. The 
even pq case is not difficult because (4-5) is essential for the equality. However for the odd 
pq case, we need careful treatment. 

Thus we will concentrate our attention on the odd pq case. Let us consider the quantity 

A{n,p, q)/A{n,p, q) and show that it is 
Further as for odd pq case, we have 



exp \ m \ l^ip + q + 1) + ^ ] {2n + qf 



1 / 1 N 

exp [ [ -{p + q + i-) + -r 



(4-8) 



/ pq p ~\~ q ~\- \ 

by considering possible cases. Noting (4-6) and (4-8), the terms exp ( ttM ^ H ^ 



multiplying with exp ( iri 
comes 



p — 1 
exp I ni 



and exp ( ni 



1 



2 2 4 



in (2-15) and (3-10) be- 
(4-9) 



As it is well-known, the Jacobi's quadratic reciprocity is given by (Ireland &: Rosen 1990, 
Chap.5) 

P\ f (l\ _ f,^ fip-'^) (9-1) 
qj \p 



exp in 



(4-10) 



It means that (4-9) expresses the Jacobi's quadratic reciprocity. Hence A{n,p, q)/A{n,p, q) 
is In other words agreement in (4-2) for the odd pq case is primitively proved. 

We note that the agreement in (4-2) essentially comes from (4-3) and (4-10) and we 
emphasize that our correspondence including the sign of the phase of Gauss sums. It is 
known that determination of the sign of the phase in Gauss sums is a very subtle problem 
(Ireland & Rosen 1990, p. 73). 

In truth, we have a natural correspondence including n/2 phase in (4-2). The phase 7r/2 
reminds us of the phase anomaly (Born & Wolf 1975 8.8.4) and caustic problem (Berry & 



Ustill 1980). (In the paper of Berry and Ustill, we can find a beauty of another connection 
between classical optics and modern mathematics.) The 7r/2 phase shift was studied in 
context of the partial differential equations as Maslov connection (Maslov 1972), number 
theory as metaplectic representations (Weil 1964), phase anomaly in optics, connection 
formula in semi-classical method in quantum mechanics, and so on. The "phase" of the 
Gauss sums is a typical example of metaplectic system. In these systems, the phase is 
a delicate object. Since the phase usually appears in optics as a phase anomaly, it is 
expected that our phase also appears due to the effects. Further it is known that in path 
integral approach, (or in this case, the Fresnel integral), the treatment of the phase needs 
more carefulness (Schulman 1981 Chap. 17; Berry & Ustill 1980; Matsutani 1997). Hence 
the agreement between both approaches from Helmholtz equation and Fresnel integral is 
highly non-trivial. 

However conversely, as we solve the same system by means of different methods, these 
physical investigations require the complete agreement between A{n,p,q) and A{n,p,q). 
It means that our computation gives a novel proof of the reciprocity of the Gauss sums 
from physical point of view. 

Even though we deal only with the classical optics, we know a following correspondence 
between classical optics and quantum mechanics. The wave length A is translated into 
Planck constant h, the Helmholtz equation corresponds to the Schrodinger equation and 
the Fresnel integral is related to path integral (Guillemin & Sterngerg 1984 1.12; Sanchez & 
Wolf 1985). In paraxial theory of geometrical optics, which is sometimes called Gaussian 
optics (Born & Wolf 1975, 4.4), we can find the symplectic structure in the angle and 
position of a ray (Guillemin Sz Sterngerg 1984 Ghap.I; Sanchez & Wolf 1985). We can 
define a Poisson bracket for the angle and the position as in the classical mechanics. In the 
sense, the wave optics is related to the quantum mechanics. In truth Hannay and Berry 
(1980) studied the Gauss sums in the quantum mechanical context, which play the same 
roles in Talbot effects in classical optics as mentioned above. 

As mentioned in the introduction, we should regard the result from Helmholtz equa- 
tion as a wavy property and one from the Fresnel integral as a particle property. The 
complementary principle leads us that complementary elements, such as configuration 
and momentum, are of dual and complementary to each other due to Planck constant h 
(Bohr 1928). Though the concept was introduced in order to explain quantum mechanical 
experiment (Bohr 1928), this concept might be applicable to an even purely theoretical 
phenomena e.g., one in number theory. In fact, in the wave expression (2-7), we find the 

term —2TTi( ^2""^^ which is proportional to A because ( — zX/a"^, whereas in the particle 

expression (3-5), ni i— — h d— + —J is proportional to 1/A. When A vanishes, former 

one vanishes while the later case diverges. (In the later case, the operation is related to an- 
other connection between number theory and quantum system (Matsutani 2001a, 2001b).) 
They show the duality and complementarity. It implies that the agreement in (4-2) can 
be interpreted as wave-particle complementarity. Thus the complementarity in this system 
and reciprocity of Gauss sums should be regarded as double aspects of a thing. 
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The complementarity between wave and particle, for the quantum mechanics, is based 
upon the commutation relation. 



xp — px — \/—lh, (4-11) 

for position and momentum operators, x and p. (4-11) is resembles to (4-3), i.e., for 

- I and I — 
<l)p I P 

-] 5-1 — 1 P=^- (4-12) 



coprime integers p and q, there exists integers { - \ and | — ^ such that 



(l)p I P ^ 

(4-12) is very essential and is the most important origin which generates a beauty of 
number theory. 

Hence this analogy between (4-11) and (4-12) leads to our conclusion that the wave- 
particle complementarity (in this system) plays the same role as the reciprocity of coprime 
numbers in number theory. 

In the studies of mathematical physics, we sometimes encounter the cases that we 
feel the resemblances between number theory and physics as mentioned in the papers 
(Matsutani 2001a, 2001b). The works of Hannay and Berry (1980) and Berry and Klein 
(1996) are the most typical cases. Further in the optical design, Gaussian optics in lens 
system can be expressed by a generalized bracket of the Gaussian bracket which represents 
continuous fraction (Tanaka 1986). In one-dimcnsional classical point-particle-systcm, the 
cyclomonic polynomial appears and expresses a sort of integrable condition (Ishiwata, 
Matsutani & Onishi 1997). Further in quantum mechanical problems, we can find the 
several resemblances. We hope that this report might have some effects on these studies 
in future. 

Appendix 

In this appendix, we will base upon the Appendix in the paper of Hannay and Berry 
(1980) and supply their arguments. For (2-15) and (3-10), we will show the derivations of 
(A-14) and (A-18), which slightly differ from those in the paper. The difference assures 
the equality in (4-2). 

First we will define an infinite sum, 

1 fin \ 

G(a, 6, c) := ^lirn^ 2^ E e^p i- [am' + cm] \ , (A-1) 

m=-Nab ^ ^ 

and the Gauss sum, 

K{a, b, c) := - ^ exp ( y [am' + cm] j . (A-2) 

^ m=o ^ ^ 
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By letting m = bn + s, each term in (A-1) becomes 



exp ( — (am^ + cm) j = exp ( — [as^ + cs) ) , 



(A-3) 



if ah and c are even (or ah and c are odd). Hence if ah and c are even (or ah and c are 
odd), (A-1) and (A-2) coincide with, 

G(a,6,c) = X(a,6,c). (A-4) 

For the case of that a and c are even and a and h are coprime, we can compute K{a, h, c) 



as. 



6-1 



ii:(a, c) = - ^ exp 

m=( 

1 

4E 



1 + 



27rza/2 
an/ 2 



c 

m+- 



exp 



—ma 



n=0 

_ 1 /^«/2 



27r?; 



—ma 



exp — — n exp 

/ \ 





2 


"1" 


.2. 




a 


1 ^ 


"1" 


D 




a 





exp 



-7rz(a/2) 



2J 



-,2\ b-1 



E 

n=0 



Jexpf— n 



(A-5) 



In the case for odd h (Ireland & Rosen 1990, Chap. 6), the Gauss sum is expressed by, 

a-1 



E 

n=0 



Jexpf— n 



for 6=1 mod 4, 
i\/h, for h — 2) mod 4. 



= v^exp(^(6-l)2 



(A-6) 



Using the relations. 



we have 



5Z (6)^^p(t^''1 = V5exp(4^(5-i; 



^ , a-1 



n=0 



27ri 



(A-8) 



Next we will consider the case that ah and c are odd and a and h are coprime. If we 
deal with exp(27rim/6)-type functions, we can use a technique mentioned in the paper of 
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Hannay and Berry (1980), i.e., X^J^^i exp(27rzam/6) = ^^pI^ti"*"^/??) for coprime 

and odd ab. However we are treating exp(7rim/&)-type functfons. Since the period of 
exp{'Kim./b) in m is not b but 26, we must independently consider this case, which is 
completely different from the even a case. 

Since 6 + 1 is even and 6 + 1 = 1 modulo 6, we have 



l + b 



(A-9) 



Hence we have 



and 



Noting (-1)' 



exp (^7ri-rnj 







"1" 


= exp ( 


27Ti 


2_ 






"1" 


= exp 




2_ 



a 2 ■ ' 
-m — TTiam' 

b 



-m — TTicm 



(A-10) 



(A-11) 



exp (- 



-mam 



+ nicrnj = 1. 



(A-12) 



Accordingly for odd ab and odd c case, we have. 



K{a, b,c) = ^Y^ exp '^^^^ 



m=0 

_ 1 fa/2 



exp 



—2'Kia 



m + c 



1 

2a 
1 

2a 

2nia, 



exp 



J 6J 

-,2\ &-1 
b/ n=0 



— 27rm 2 1 
2 

n\ [2111 



2a 



1" 




" 1 " 




2 


6 


_2a_ 


:) 



(A-13) 



Here following the argument in the paper of Hannay and Berry (1980), we will prove 
the reciprocity of the infinite sums for any a, b and c. 



na 



G(a, b, c) = G(—b, a,c)\ — exp 





Aab 



(A-14) 



We first note the relations; 



^ J 

12 



/3 



(A-15) 



CO 




= J2 e^'^*^^. (A-16) 

m= — co 

Secondly, we note that due to the paper of Hannay and Berry (1980), the left hand side 
in (A-14) is expressed by, 

G{a, b, c) = lim [ v^e-'^^^'^clx) exp ( ^(aa;^ + ex)] dx. (A-17) 

Using (A-15) and (A-16), we compute (A-17) and then derive (A-14). Using (A-4), (A-14) 
becomes the reciprocity of the Gauss sums. 
Hence the Gauss sums are expressed by. 




a even, h odd, 
h odd, a even, 
a odd, h odd. 



(A-18) 
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